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Abstract 

Graviational radiation is described by canonical Yang- 
Mills wave equations on the curved space-time mani- 
fold, together with evolution equations for the metric 
in the tangent bundle. The initial data problem is 
described in Yang-Mills scalar and vector potentials, 
resulting in Lie-constraints in addition to the familiar 
Gauss-Codacci relations. 

1 Introduction and outline of 
the approach 

The asymptotic gravitational wave structure resulting 
from the coalescence of astrophysical black holes or 
neutron stars is receiving wide attention in connec- 
tion with the gravitational wave detectors currently 
under construction. A compact binary system pro- 
duces strongly nonlinear gravitational waves of which 
a small residue escapes to infinity as radiation. The 
asymptotic wave structure at a distant observer is in 
quantitative relation to the system parameters of com- 
pact binaries, which makes gravitational radiation a 
new spectrum for astronomical observations. The pre- 
diction of gravitational wave forms is presently pursued 
by numerical simulation. The current approaches are 
based the 3+1 Hamiltonian formulation by Arnowitt, 
Deser and Misner [^J, with the notable exception of 
strategies employing null-coordinates (e.g. 0) or sys- 
tems of conservation laws pj . 

The significance of gravity waves in compact bina- 
ries, both in the process of coalescence and in radia- 
tion, suggests to focus on a description of relativity 
by nonlinear wave equations. Wave equations pro- 
vide a proper frame- work for the physics of nonlinear 
wave motion, establish connection with electromagnet- 



ics (abelian and linear) and Yang-Mills theory (non- 
abelian and nonlinear; see, e.g., and may comple- 
ment present approaches in numerical relativity by way 
of understanding or circumventing some of the compu- 
tational difficulties. This paper establishes a first step 
in this direction. Generally, wave equations provide a 
setting appropriate for numerical implementation via 
cylindrical or spherical coordinate systems, including 
outgoing boundary conditions at grid boundaries, and 
offer a possible setting for numerical treatment of in- 
going horizon boundary conditions. 

Pirani Q presented compelling arguments to con- 
centrate on the Riemann tensor in describing gravita- 
tional radiation. One can proceed by proposing to es- 
tablish a numerical algorithm based on the divergence 
of the Riemann tensor R a bcd, in which the tensor 



Ucd = ^(VfcT^f, - -gb[d^c]T) 



(1) 



acts as a source term. A tetrad {(e M ) ft }-formalism 
takes the Riemannian formulation into equations of 
Yang-Mills type. With Lorentz gauge on the connec- 
tion 1-forms, uj n 







(2) 



(V Q denotes the covariant derivative associated with 
the metric, g a b), fully nonlinear, canonical Yang-Mills 
wave equations are obtained 



(3) 



Here, □ is the non-abelian generalization of the 
Laplace-Beltrami wave operator. The wave equations 
are also derived for the scalar, Ricci rotation coeffi- 
cients. Complemented by the equations of structure 
for the evolution of the tetrad legs, a complete system 
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of evolution equations is obtained which is amenable 
to numerical implementation. 

This interwoveness of wave motion and causal struc- 
ture distinguishes gravity from the other field theories. 
In the present description, this two-fold nature of grav- 
ity has thus been made explicit. Essentially, gravita- 
tional waves are now propagated by wave equations on 
the (curved) manifold, while the metric is evolved in 
the (flat) tangent bundle by the equations of structure. 
Of course, such two-fold description is only meaningful 
for wave motion with wave lengths above the Planck 
scale, below which the causal structure is not well- 
defined and wave motion can not be distinguished from 
quantum fluctuations. 

Because the wave equations are second order in time, 
conditions on the initial data involve both U) afil , and 
their normal Lie derivates L n L0 a(LV . In addition to the 
Gauss-Codacci relations, w aAI „ and C n Lo a ^ v are related 
to the initial energy-momentum and Tb^v The latter 
follow from Gauss-Riemann relations, as non-abelian 
generalizations of Gauss's law. 

2 Equations for R a bcd 

We will work on a four-dimensional manifold, M , with 
hyperbolic metric g a b- In a given coordinate system 



{x } the line-element is given by 
ds 2 = g ab dx a dx b . 



(4) 



The natural volume element on M is e a tcd = 
y/—g~A. a bcdi where g denotes the determinant of the 
metric in the given coordinate system, and A a bcd de- 
notes the completely antisymmetric symbol. Covariant 
differentiation associated with g a b will be referred to by 
V a . The geometry of M is contained in the Riemann 
tensor, R a bcd, which satisfies the Bianchi identity 



3V[ e -R a 6]cd — ^ c Rabcd + ^ aRbacd 
+VbReacd = 0. 



(5) 



Using the volume element e a bcd, the dual *R is defined 



as 



■ef 



form 



2 ab 



R 



efed- 



The Bianchi identity then takes the 



V a * Rabcd = 0. 



(0) 



Einstein proposed that the Ricci tensor R a b = R t 



acb 



and scalar curvature R = i? c c are connected to energy- 
momentum, T a b, through 



1 

Rab — 1^9abR = 87rT a fc. 



(7) 



Geometry now becomes dynamical with (0) (see, e.g., 

H) 



V Rabcd = 2V[f,-R a ], 



(8) 



In the presence of Einsteins equations (0), therefore, 
(||) becomes 

E : V a Rabcd = 16TT(V [c T d]b - \g b [d^c\T). (9) 

The quantity on the right-hand side shall be referred to 
as Tbcd- In vacuo, E has been discussed by Klainerman 
S, who refers to E (with Tbcd = 0) together with (||) 
as the spin-2 equations. The tensor Tbcd is divergence 
free: 



V b T bc d = 0, 



(10) 



in consequence of the conservation laws V a T a {, = 0, 
and in agreement with V h V a R a b c d = 0. 

2.1 Yang-Mills equations 

In the language of tetrads, additional invariance arises 
due to the liberty of choosing the tetrad position at 
each space-time point. This invariance is described 
by the Lorentz group, and introduces a vector gauge 
invariance in the form of the connection 1-forms, gov- 
erned by Yang-Mills equations as outlined below. The 
formal arguments can be found in the theory of super- 
symmetry (see, e.g., |], |, §). 

Let {(e^) 6 } denote a tetrad satisfying 



riiiu = <^sign(77 Ml/ ). 



(11) 



Neighboring tetrads introduce their connection 1- 
forms, w apv , 



(e^Vaie^c 



(12) 



The connection 1-forms uj aiiv serve as Yang-Mills con- 
nections in the gauge covariant derivative 



(13) 



satisfying V a (e M ) b = 0. Here, the commutator is de- 



V a = V a 
0. 

fined by its action on tensors 4> ai ---a k a 1 ---a l 

l^ni 4 > a 1 ---a k ]a 1 ---cti 

ECtj 
i U>aa i LOa 1 ---a k ai---a j ---ai- 



as 



(14) 
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In particular, we have In (p3|), 9t(e Al )t is left undefined. Denning £ h = (dt) b , 

the four time-components 

[L) a , Wfejjui/ = UJ al ^UJbaL> — ^> a v' L0 ba i i- (15) 

In what follows, Greek indices stand for contractions 
with tetrad elements: if v b is a vector field, then = 
v b (e fl ) b , and = ff v v v 



:= (e^) a C (24) 
become freely specifyable functions. The evolution 



fj_)b, anu v i\ u v . equations for the tetrad legs thus become 

The Yang-Mills construction thus obtains for the 



dtie^t, + io»(e„) b = d b N^ + u> b "N v . (25) 



Bianchi identity 

abf * u ' ^ ^ The tetrad lapse functions iV^ are related to the famil- 

and for the equivalent of E mr ^ a P se an d shift functions in the Hamiltonian for 



E' : V a R a btiu = Thuv (17) 



malism through 

gat = N a (e a ) a = (N q Ni - N 2 ,N p ). (26) 



The Bianchi identity ( |l6| ) introduces the representation 

(cf. Jl3)) We will now turn to evolution equations for the con- 



Rabp.v = V 1 aLOh^v ~ ^b^a^u + [^a,^b\nv, (18) 



nection 1-forms. 



which will be central in our discussion. o "c 1 j. • c 

The antisymmetries in the Riemann tensor introduce HjCJUatlOIlS IOF (jJaiiv 

conditions on initial data on an initial hypersurface, S. 

If v b denotes the normal to £ and ^ e define a Lorentzian cross-section of the tangent 

bundle of the space-time manifold by [[ll] 

:= V d u d ^ = 0. (27) 



V„=^„KV c )+D a onS, (19) 
we obtain the Gauss-Riemann relations 



V b D a R a bcd — ~Pcd 
,b na 



The Lorentz gauge (27) provides^] a complete, six-fold 

, ( (20) connection between neighboring tetrads. The six con- 

v * abed ■ straints = are incorporated in E' by application 

, h _, , . , . of the divergence technique JTol, O : 

where p c d = v T hcd . Inese may be considered gener- H — 1 1 — H 

alizations of Gauss's law in electromagnetism. Condi- w i, fjarn , i_ 

tions © find their equivalents in the tetrad formula- : 1 ab ^ + 9abC ^l ~ Tb ^ W 

tion. To this end, we write, analogous to ([l9|), on £ 

the derivative as Recall the transformation rule for the connection 

(see, e.g. @) 

V a = -v a {v c V c )+b a . (21) 

w aP B = A^AjAw + AfVaApa. (29) 

The Gauss-Riemann relations (EOI) thus become 



In the present tetrad language, this gauge trans- 
v D a R a b^u = p^v, ,„_.. formation is readily established by consideration of 

v b t) a * Rab^u =0. two tetrads, {(e M ) b } and {{e f _ l ) b }. The construction 

Ap v := (ep) c (e ly ) c provides a finite transformation 
2.2 Evolution equations for the tetrads v fi = ^ a v a of a field v a = (e a ) b v b in the {(e M ) b }- 

tetrad representation into v a in the {(e M ) 6 }-tetrad rep- 
The definition of theconnection 1-forms gives the equa- resentation This applied to pj obtains (Re* _ 



tions of structure lie 



1 Conceptually, c M „ = f(i^ a ^v,9ab) w iU a l so serve its purpose, 
^[a(^/j)i>] = ( e!y )[fc w a]i/^i- (23) where /(-,-) depends analytically on its arguments. 
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3.1 Existence of Lorentzian cross- 
section 

To proceed, we consider in an open neighborhood 
A/"(E) of E with Gaussian normal coordinates {r, x p } 
(i> c V c t = 1) the infinitesimal Lorentz transformation 

A/ = 5/ + ir 2 <7/m^(E). (30) 



It is convenient to employ language of scalar and vector 
potentials, and A a ^ Vl respectively, defined in 

W aful = Va^^y + A afl „. (31) 

Denoting the effect of (|o|) via ( p9| ) by a superscript 
(r), we have 

= ^ + r<S> MJ/ in Af (E), (32) 



so that 



=^-^ inAT(E). (33) 



Using a geodesic extension of the normal ^ b off E, con- 
straints (p7j) become 



In A„ 



pv = VcUJ^ = + D c w. c MV on E. (34) 
this obtains 



(35) 



<V = "V + D c (v c $ MV + A c ^) 
= $ M „ + + DoAf^ 

on E. Consequently, we have the transformations 

r"~ I 

on E. (36) 



,4 (r) - 4 

By choice of cr^^ = c Mt , , the constraints 
into 



0. 



transform 



(37) 



Notice that bringing the tetrad in Lorentz gauge (|3J 
by ([30|) is achieved by proper second time-derivative 
d^{E^ T=0 of its legs, leaving the tetrad position and 
its first time-derivative as invariants at r = 0. 

We shall now establish that E" maintains the 
Lorentz gauge ( [27] ) in the future domain of depen- 
dence of E. The inhomogeneous Gauss- Riemann re- 
lation (^0|) is implied by antisymmetry of the Riemann 
tensor in its coordinate indices, and gives 

= V b {V a (R abf _ lL , + QabC^u) - n^u} 

= v h {D a R ahtlv - T htiv ) + (i/ 6 V 6 )c p „ (38) 
= (^V b )c^. 



The inhomogeneous Gauss-Riemann relations are 
gauge covariant, so that we are at liberty to consider 
initial data satisfying both (|3^) and the gauge choice 
d37T), whence 



c^v = (^ c V c )c M „ = on E. 



(39) 



In (|28|) c p „ satisfies a homogeneous Yang-Mills wave 
equation (c/. 0, HI), 



ticuv ~ V c V c cw = 0. 



(40) 



It follows that the scalar fields c^ v are solutions to an 
initial value problem for homogeneous wave equations 
(flc|) with trivial Cauchy data (^9|). Consequently, 



= in D+(E), 



(41) 



which establishes that solutions to E" are solutions to 
E'. 

We now elaborate further on (p8|). 
3.2 Canonical wave equations 



In Lorentz gauge 02 7|) , the divergence equation E" 
obtains wave equations for the connection 1-forms 
through the representation of the Riemann tensor (|l8|). 
Indeed, by explicit calculation, we have 



(42) 



Here, we have used c M „ = 0, so that V 'a c nv = V a c MI/ . □ 
is used to denote the Yang-Mills wave operator V C V C . 
The Ricci tensor R a b in ([l2]) is understood in terms of 
T a b using Einstein's equations. 

Similarly, we can obtain a system of scalar equations 
for the Ricci rotation coefficients, Lo a ^ v = {e a ) a uj a ^ v . 
This involves a number of manipulations, the result of 
which is 



(43) 



The two-fold nature of gravity can now be expressed 



as 



Separation Theorem. 1 Gravitational wave motion 
is governed by canonical wave equations on M . In re- 
sponse to the wave motion, the metric on M evolves in 
the tangent bundle of M by the equations of structure. 
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The coupling of matter to the connections is ac- 
counted for by Ttcd- It is of interest to note that T bc d 
contains vorticity; a detailed enumeration of the nature 
of Tbcd falls outside the scope of this dicsussion. 

The Theorem suggests some computational approx- 
imations for weakly nonlinear gravity waves. Firstly, 
consider equations with uncoupled (prescribed or 
fixed) metric, with given Laplace-Beltrami wave op- 
erator, D g . Thus, we have (in vacuo) 



S W a/1 „ - [LU C , Va^c]^ = 0, 



(44) 



with the metric g a b as a background field providing the 
underlying causal structure. The approximate radia- 
tion then follows from the fluctuations in the metric as 
would follow from the equations of structure. Secondly, 
small amplitude waves allow us to linearize, thereby 
obtaining the purely abelian wave equations 



^riUauv — 0, 



(45) 



where rj refers to the Lorentz metric. Numerical exper- 
iments must show the validity of such approximations. 

4 Initial value problem 

Initial data for the wave equations are L0 afll , and its 
Lie derivative £ n Lu afJ-l , on E. These data must satisfy 
certain constraints on E, commensurate with the ini- 
tial distribution of energy-momentum and the Gauss- 
Riemann equations (22). We shall express these equa- 
tions in terms of scalar and vector potentials (|l]) . 
The projection tensor, h ab , onto S is 



g a b + v a v b . 



(46) 



The covariant derivative induced by h ab shall be de- 
noted by D a , i.e., D a h c d = 0. In what follows, 
(i/ c V c )/ is also denoted by /. The unit normal is ex- 
tended geodesically off E, so that v c V c v b — and 
(v c V c )h ab = h ab {v c V c ). Thus, $ = (^ C V C )$ and 
A afiV = {v c V c )A allv . We further define 



A'. 



K*, 



— r A —A i ft- M v ^ 



4.1 Expressions for A 



afiu 



Consider a normal tetrad, {(£' /i ) f> }, in which one of 
the legs is (initially) everywhere normal to the initial 
hypersurface: 



v b on E. 



(48) 



Extrinsic curvature, K ab , of E can be defined 'static' by 
projection of the unit normal following parallel trans- 
port over E, 



K a b = D a v b , 



(49) 



and 'dynamic' by the Lie derivative, \L n h abl of the 
projection operator h ab with respect to a Gaussian nor- 
mal coordinate, n. With a normal tetrad, its (extrin- 
sic) helicity, i.e, the twist in a strip swept out by the 
integral curves of a leg (E u ) b passing through an inte- 
gral curve of leg (£^) b is 



En 



Vb(E ll ) c V c {E v f on E. 



(50) 



For /x, v ^ n, describes twist in E when fj, ^ z/, 
and bending when fj, = v. It follows that 



H a0 h a X = -K, 



(51) 



in view of (|48j): v a (E a ) a = -8% on E. The normal 
(fi = n) and tangent (/i ^ n) extrinsic helicities are 



H^v = u b (E^) c {uj CT (E'') b } 
-<5> nu if /i = n, 



(52) 



By (51), the symmetry of K ab and A^ nl , 
have for fj,,u ^ n the symmetry 

A[ij,v] n = 0. 

If fj, n, (E^a = h ab (E fJi ) b , and so 

A aiiv = h b a (E ll ) c V b (E„) c 

= {E ll ) c [h d h b a V h ]{E v y 
= (E^) c D a {E u ) c . 



-Aftvn, we 



(53) 



(54) 



If also v 7^ n, D a acts in (54) on tangent legs (E v ) b , 
the result of which is determined by h ab . We shall 
write A aiiv — A aa ph^h^ for this intrinsic part of A allv . 
In the normal tetrad, therefore, A aiiv falls into two 
groups: (i) the extrinsic part A^ vn — K^, and the (ii) 
intrinsic part A ajLV1 with which we have 



+ 2K a[ ^ } on E. 
4.2 Expressions for R a b^u 



(55) 



On E, R a bu.v contains intrinsic, three-dimensional, and 
extrinsic parts by substitution of (|3l|): 



Rab^u — V a UJ bf j, u — VfeOJa^i/ + [ui a , CU b ]iJ,v 

= 2v [b b a] ^^ + 2v [b A a] ^ (56) 
+2D[ a Ab] MV + [AajAblfiv. 
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Using (p5|), the latter may be reduced to 



By Einsteins equations, we have 



Rabfiv — 2!/[ 6 {£)a]$ Ml/ + A a ]iJ,u} 

+2D [a A b]liV + [Aa,At]p, (57) 
+4D [a K b][ ^ ] + 2K a{fl K u]b , 

where b a K^ b = D a K^ b + A a ^K lb . Notice that the 

intrinsic, three-curvature ^' R abc d of £ associated with 
h ab satisfies 

{3) Rab^ = D a A hiy - T) b A ailv + [A a ,A b ]^. (58) 
It follows that 



Kftl Rcdfiv — 2D[ a Aqn V + [A a ,A b ] 



+2K a[ ^K ly]b + 4D [a K b]l ^ ] 
( ^R ab ^ + 'ib [a K b][ ^ ] 



(59) 



This last equation generalizes a similar expression ob- 
tained in j|, where fi,u ^ n is considered. 

4.3 Energy-momentum relations 

Constraints related to the initial distribution of energy- 
momentum are obtained by consideration of the Ricci 
tensor. To this end, consider the two expressions (|5£" 
and ( |59| ) for Rab^v of the previous section in 



h c a V RcbfiV — DaQ^u A' , 



hahfrRcdfm = 2D[ a K b ] ll . 



(60) 



Consequently, the elements of the Ricci tensor become 



' J R bu = v b h cf *R, 



h d b R dn = 2h d b h a »D [a K d]tl 
= D a K ab -D b K. 



(61) 



The second expression (q9|) for the Riemann tensor also 
also gives 



h«hflh C a hiRa ba = (3) Rabcd + K atl K vb 
KaisR [ib: 

resulting in the familiar identity 

167rT„„ = 2{R ab ~ \g ab R)v a v h 

— tl tl Kabcd 

= h a »h b PR abaf3 

= WR + K 2 - K ail K a ^ 

= ^R + K 2 - K ab K ab . 



(62) 



(63) 



R ab = Sn(T ab - ^g ab T) = 8wf ab . (64) 

Combining the results above, we have the constraints 
on the Lie derivative of A ativ in 



(65) 



together with the familiar Gauss-Codacci relations 



WR + K 2 - K ab K ab = 167rT„„, 
D a K ab -D b K =8nh%f an . 



(66) 



We shall refer to (|65|)-(|66|) as the energy- momentum 
constraints. 

In going from a first order Hamiltonian description 
to a second order description, it is not surprising to en- 
counter also constraints on the Lie derivative, reflecting 
the conditions that the Gauss-Codacci relations also 
must be satisfied on a future hypersurface in 

4.4 Gauss-Riemann relations 

The Gauss-Riemann relations (^2|) can similarly be ob- 
tained in terms of (fr^ and A aliv . The antisymmetry 
of the Riemann tensor in its coordinate indices implies 

v b D a R abfJil/ = D a (y R ablJiV ) — R ablJjV K a , . 

= D-(y b R abtiV ). [b '> 



By (|56D, it follows that 

V Rab^LV ~= Da^tiv Aa^ J^.^ A c ^ v 

— A — K A 

= —Daunts — A' aflly . 



(68) 



Therefore, the first, inhomogeneous Gauss-Riemann 
constraint in (E2) reads 



v b D a R abilv — v b D a R abil _ v + [A a , v b R ab \ jlv 
= D a {v b R abllv ) + [A a ,v b R ab ] llv 
= b a {v b R ab ^) 



(69) 



where the symmetry of the extrinsic curvature tensor 
was used in going from the third to the fourth equation. 
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These six constraints on C n A a ^ u arise in view of the 
functional relationship 



(70) 



Because \C n h ab = K ab and K afl = A afin , six degrees 
of freedom in 



(71) 



are constraint by (JTOj) . 

To summarize, the constraints on the initial data in 
terms of potentials are 



Gauss- Codacci: 



WR + K 2 -K ab K ab = 16ttT„„ 
D a K ab -D b K = 8nh^f an 



Lie- constraints 



iiv -r <™-nv — —8nT v bi> 



A% 

D a D a <S>^ + D a A' aflu = - P|t „ 
Lorentz gauge: 



+ D a A afW = 0. 



(72) 



(73) 



(74) 



The condition on a'" ensures that the second Lie- 
constraint remains satisfied. This shows that by choos- 
ing 



a,,„ = -K<b, 



we obtain 



= on E. 



(78) 



(79) 



With this restricted Lorentz gauge, the Lie-constraints 
reduce to 



D a A 

+ b c A° 



ativ 



0. 



on E. 



(80) 



We remark it does not appear to be feasible to en- 
sure $ = throughout D + (E) by suitable choice of 
restricted Lorentz gauge on the initial data. 
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4.5 Restricted Lorentz gauge 

Restricted gauge transformations for the connection in 
analogy to the same problem in electromagnetics, can 
be considered through the Lorentz transformation 



A; = + tA; + in Af (E), (75) 

where r is a Gaussian normal coordinate of E such that 
v c V c t = 1, and d T X^ = 0. By @, we have 



A {r) - A 



> on £. 



(76) 



The first (energy-momentum) Lie-constraint is gauge 
invariant (invariant under transformations ( |29| ) ) , while 
the second is not. Therefore, we can choose X" and 
<j,Y so that 



4 r i = v-v = ° 



(<£ (r) )U 



on E. 



(77) 
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